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Abstract 

The energy loss of a hard parton is regarded as a useful tool to probe properties of 
a soft medium. However, the information obtained by a measurement of the nuclear 
suppression Raa is not very sensitive to medium properties. Thus, more differential 
observables are needed. Here we focus on the pattern of angular correlation of 
hadrons associated with a hard trigger. These show rich pattern when going from 
low pt to high pt associate momentum. At low px, the pattern of associate hadrons 
can be understood in terms of Mach shocks excited by the energy lost from the away 
side parton while the parton itself is absorbed by the medium. At high px, fluid 
dynamics is not applicable any more and the correlations become dominated by the 
punchtrough of the away side parton with subsequent fragmentation. We present 
an analysis of these phenomena taking into account the full expansion dynamics of 
soft matter and show how both transverse and longitudinal expansion play a crucial 
role in understanding the data. 
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1 Introduction 



The expression 'jet tomography' often used to describe the analysis of hard 
pQCD processes taking place inside the soft matter created in an ultrarela- 
tivistic heavy-ion collision. In particular the focus is on the nuclear suppression 
of hard hadrons in A-A collisions compared with the scaled expectation from 
p-p collisions due to loss of energy from the hard parton by interactions with 
the soft medium (see e.g. [1]). However, the nuclear suppression factor 

d 2 N AA /dp T dy 

Raa{Pt > V) ~ T AA (h)d^ NN /d PT dy- (1) 
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is a rather integral quantity, arising in model calculations from the schemat- 
ical convolution of the hard pQCD vacuum cross section da^.~^^ +x for the 
production of a parton /, the energy loss probability Pf(AE) given the ver- 
tex position and path through the medium and the vacuum fragmentation 
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Here, fi/A(%,Q 2 ) denotes the nuclear parton distribution function dependent 
on the parton momentum fraction x and the hard momentum scale Q 2 and 
f + k is the the partonic pQCD cross section. 
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Eq. (2) has to be properly averaged over all possible vertices distributed ac- 
cording to the nuclear overlap Taa and paths through the medium. In [2] we 
have argued that it is possible to factorize this spatial averaging from the 
momentum space formulation Eq. (2) and thus define the geometry-averaged 
energy loss probability (P(AE, E))t aa - Raa can thus be viewed as providing 
constraints for the form of (P(AE, E))t aa - 
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Fig. 1. Left panel: Geometry-averaged energy loss probability distributions for differ- 
ent energy loss scenarios [2]. Right: Nuclear suppression factor Raa for the different 
scenarios compared with PHENIX data [3]. 

In Fig. 1 we show trial distributions (P(AE, E^}t aa based on different possible 
scenarios how the medium induces jet quenching. Although the form of the in- 
dividual (P(AE, E)) Taa differs substantially, the resulting Raa is very similar. 
This implies that Raa does not provide strong enough constraints to extract 
medium properties reliably except for the fact that quenching is substantial 
[2]. However, given that assumptions about the longitudinal and transverse 
flow profile (while requiring that the soft hadron spectra agree with data) can 
alter the medium quenching power by a factor of five [4] , this means that any 
attempt to gain access to information about the medium density must involve 
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observables beyond Raa- In [2] we demonstrated that 7-hadron correlations 
[5] are capable of resolving the differences; in this paper we're interested in 
the additional information accessible by angular correlation measurements. 



2 The framework 

Energy and momentum lost by a hard parton reappear in the medium. We 
make the assumption that the bulk of the lost energy excites a shockwave in the 
soft medium. If we study back-to-back production of hard partons and require 
one of the partons to yield a hard hadron trigger, the angular correlation 
pattern associated with this trigger consists of two components with distinct 
properties. First, there are the remnants of the hard processes — next to 
leading hadrons created in the fragmentation alongside the trigger and hadrons 
created in the fragmentation of the away side parton whenever this traverses 
the medium. The only place medium properties appear in this contribution 
is the amount of energy lost from the partons, but as long as fragmentation 
occurs in vacuum, neither correlation angular position nor width is determined 
by the medium, only the correlation strength. In contrast, of the soft process 
(i.e. the excited shockwave) only the strength is not determined by the medium 
but by the lost energy. The position is largely determined by the medium speed 
of sound c s [7,8] and its width by the freeze-out conditions [8]. 

We describe the hard underlying process by Monte-Carlo sampling Eq. (2) for 
the near and away side. The energy loss probability P(AE) given the path is 
determined by using BDMPS radiative energy loss in the formulation of [9] 
by evaluating line integrals over the local quenching power along the path of 
the hard parton [4,6]. As a description of the medium we mainly consider a 
parametrized evolution capable of describing a plethora of bulk observables 
[10] but we also study a hydrodynamical evolution [11]. 

In order to simulate the shockwave we follow the flow of energy and mo- 
mentum. We assume that a fraction / of the energy and momentum lost 
to the medium excites a shockwave characterized by a dispersion relation 
E = c 2 s p and a fraction (1 — /) in essence heats the medium and leads to 
collective drift along the jet axis to conserve momentum [7]. We calculate c s 
from a quasiparticle description of the equation of state. The dispersion rela- 
tion along with the energy and momentum deposition determines the initial 
angle of propagation of the shock front with the jet axis (the 'Mach angle') 
as = arccosc s . We discretize the time into small intervals At, calculate 
the energy deposited in that time as E(t) = At ■ dE/dT and then prop- 
agate the part of the shockfront remaining in the midrapidity slice (i.e. in 
the detector acceptance). Each piece of the front is propagated with the lo- 
cal speed of sound and the angle of propagation is constantly adjusted as 
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= arccos JJ c s (t)cIt/ (t — te). Once an element of the wavefront reaches the 
freeze-out condition T = Tp, a hydrodynamical mode cannot propagate fur- 
ther. In the local restframe, we have a matching condition for the dispersion 
relations E = c 2 s p and E = y/W+p 2 - M where M = V (p(T F ) + e(T F )) is 
the 'mass' of a volume element at Tp. Once we have calculated the additional 
boost <u* hocA: a volume element receives from the shockwave using the matching 
conditions, we use the Cooper-Frye formula to convert the fluid element into 
a distribution of hadrons. 



3 Low px associate — recoil of the hydro medium 

At an associate cut of about 1 GeV, the system is clearly dominated by the 
recoil of the hydro medium, we observe almost no away side partons emerging 
from the medium in the simulation. The shockwave shows a strong influence 
on longitudinal and transverse flow: Due to the fact that the rapidity y of 
the away side parton is not fixed by the rapidity of the trigger hadron, an 
averaging over y leads to an apparent shrinking of the cone [13], cf. Fig. 2. 
For a scenario which does not couple to longitudinal flow, this effect leads to 
disagreement with the data. 



Fig. 2. Calculated 2-particle correlation under the assumption that a) the away side 
parton is always at midrapidity and the excited mode doesn't couple to flow (green) 
b) using realistic P(y) and assuming that the excited mode doesn't couple to flow 
(red) and c) including realistic P(y) and longitudinal flow elongation compared with 
PHENIX data [14]. 

However, for a hydro mode the shock cone is longitudinally elongated due 
to its coupling to flow, therefore its spatial extension is the solution of the 
characteristic equation 
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and the subsequent widening in rapidity space leads to agreement with the 
data. Any scenario which does not couple to longitudinal flow must thus nec- 
essarily exhibit an even wider angle before averaging. 

p T = 1 GeV p T = 0.5 GeV 
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Fig. 3. The effect of transverse flow on the observable correlation peak for configura- 
tions in which transverse flow and cone are aligned (solid) and orthogonal (dashed, 
6 = 3 fm) for associate hadron = 1.0 GeV (left panel) and 500 MeV (right panel). 
All calculations are for fixed energy E, intrinsic kx = and y a way = 0. Correlation 
strength in the direction of the away side parton is suppressed for clarity. 

Transverse flow has no great influence on the position of the angle, but due 
to the fact that high p? correlations in a hydro + shockwave description are 
only generated when flow and shockwave are almost aligned, it may lead to 
the apparent disappearance of one wing of the shockwave at sufficiently high 
Pt (which is recovered at lower pt for some initial vertices. We illustrate this 
in Fig. 3. This has considerable impact on the interpretation of 3-particle 
correlation data. 



4 Raising — competition between hard and soft mechanisms 



It is clear that the relative importance of the hard and soft component of the 
correlation will depend on both trigger threshold and associate hadron cut. 
As soon as hydrodynamic flow of thermalized matter ceases to be the chief 
mechanism of hadron preduction at given px, it cannot longer be assumed that 
Shockwaves dominate the correlations. Even for a 20 GeV trigger hadron, one 
cannot expect to find shockwave correlations at 6 GeV, simply because there 
is no hydro medium to support them at this scale. Thus, the rising associate 
cut will in general shift the focus from soft to hard modes. The expected effect 
of a rising trigger is somewhat more complex. Increasing the parton energy 
will increase the likelyhood of observing punchthrough, but also the yield on 
the near side (caused by next to leading fragmentation) as compared to the 
away side (cause by Shockwaves) will change with the excitation function of 
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these mechanisms with energy. Hard processes are expected to dominate at 
large enough scales. 

Both the vanishing of the low p T dip (possibly due to punchthrough) and the 
drop of the away side orrelation strength compared to the near side expected 
from the model are qualitatively seen in the data [15] for increased trigger 
and associate momentum scale. Especially the fact that the away side pattern 
remains broad for increasing associate momentum and only drops in strength 
is a strong point in favour of a hydrodynamical mode, as pQCD motivated 
emission scenarios predict a shrinking of the angle. However quantitative re- 
gions in the transition regime between hard and soft physics are challenging 
and have not seriously been attempted so far. 



5 Hard trigger and hard associate — punchthrough 



If the trigger px is raised in the region above 6 GeV, punchthrough of the 
away side parton through the expanding medium is observed, and the yields 
obtained in the model agree well with the data, cf. Fig. 4, except in the 4-6 
GeV away side momentum bin. This indicates that the hadron production in 
this region in the presence of a medium is not primarily given by fragmentation 
of a hard parton but that sizeable contributions of other processes are present, 
quite possibly recombination [16] plays a role here. 
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Fig. 4. Yield per trigger on the near side (left panel) and away side (right panel) 
for trigger hadron above 8 GeV in the model calculation as compared with STAR 
data [18]. 

In Fig. 5 we show the distribution of vertices leading to a triggered hadron 
(note that this distribution is also relevant for Raa- While there is some de- 
gree of surface bias visible, the core of the system is by no means black and 
contributions to Raa can come also from the core of the matter. This is in 
contrast to previous studies [17] and the difference cen be traced back to the 
fact that the full expansion is taken into account in the present model [19]. 
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Fig. 5. Normalized density of vertices for events with a trigger hadron above 8 GeV. 
The near side parton propagates into negative x direction. 

In Fig. 6 we illustrate the spatial region probed in back-to-back dihadron 
measurements. This region shows rather large model-dependence, but is in 
general peaked more towards the medium center. This geometry dependence is 
a promising handle to gain information about the medium density distribution 
from dihadron correlation measurements [19]. 



Hydrodynamics 



Box density 





Fig. 6. Normalized density of vertices for which an associate hadron in the momen- 
tum region 4 < pt < GeV is found, given a trigger in the region 8 < pr < 15 GeV. 
The near side parton propagates into negative x direction. 



6 Summary 



We have argued that Raa is not sufficiently sensitive to medium properties 
and that more differential measurements are needed. We have presented a 
consistent picture of angular correlation pattern given a hard hadron trigger 
based on hard processes and the recoil of the soft medium. Within this unified 
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description, the correlations probe both the medium flow profile and the den- 
sity distribution, thus in principle more detailed information can be gained. 
In practice, more high-precision data are necessary — while it is easy to rule 
out toy models, the sensitivity to distinguish different dedicated bulk matter 
models is not yet sufficient [19]. 
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